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I propose a simple and manageable method that allows for deriving coupling constants of model 
energy density functionals (EDFs) directly from ab initio calculations performed for finite fermion 
systems. A proof-of-principle application allows for linking properties of finite nuclei, determined 
by using the nuclear nonlocal Gogny functional, to the coupling constants of the quasilocal Skyrme 
functional. The method does not rely on properties of infinite fermion systems but on the ab initio 
calculations in finite systems. It also allows for quantifying merits of different model EDFs in 
describing the ab initio results. 
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Approaches based on the density functional theory 
(DFT) provide us with very efficient and useful tools 
to describe properties of many-fermion systems like 
molecules, solids, or nuclei 0-0. The range of possible 
various applications and implementations is extremely 
wide. In electronic systems, there exist numerous meth¬ 
ods and techniques of linking the density functionals to 
the underlying Coulomb interaction, but in nuclei such 
links are much more difficult to explore, primarily be¬ 
cause of the fact that nucleon-nucleon interactions are 
less obviously definable. 


Only fairly recently, families of interactions based on 
chiral effective field theory have become a gold 

standard for ab initio calculations of nuclear proper¬ 
ties [l3-E3|- Based on this developments, there were 
already many attempts to link the nuclear DFT to 
these ab initio approaches 2Cj-27|. Up to now, most 
of them referred directly (or indirectly, through the so- 
called density-matrix expansion (DME) [28]) to infinite- 
nuclear-matter properties. In this paper, I propose 
a generic method that directly links the ab initio ap¬ 
proaches to the DFT in finite fermion systems. 


A classic formulation of the DFT relies on a variational 
approach to the many-body problem, whereby it assumed 
that we are able to perform an exact variation of the 
average energy 6E = <5(T|IT|'k) = 0, which gives us the 
exact ground-state energy Eq and the exact ground-state 
wave-function |To). By replacing the full variation with 
a two-stage variation, the DFT then appears in a very 
natural way. 


The first-stage variation is performed under the con¬ 
straint that the one-body local density of the system, 
p{r) = (\I , |a+a T .|\f r ), is fixed to a given density profile. 
Performing such a variation for all density profiles, one 
obtains the exact energy density functional (EDF) E \p\. 
In principle, such a constrained variation can be realized 
by an unconstrained variation performed for the system 
placed in an external one-body local potential —U(r), 


that is, 


SE' = 6('S?\H-U\'&) = 5 


('Pl./Tl'P) — drU(r)p(r) 


= 0 , 
(!) 

whereupon the external potential acquires a role of the 
Lagrange multiplier. In a sense, the constrained variation 
corresponds to probing the system with an external one- 
body field. By inverting the obtained relation p [[/] and 
inserting it into the functional E \JJ\ , one can, again in 
principle, obtain the final exact EDF E \p\. These classic 
arguments have the same structure as the effective-action 
approach, see, e.g., Refs. 21, IMS]. 

The second-stage variation, with respect to the density, 
5 p(r) E \p\ = 0, obviously then gives the exact ground- 
state energy E 0 and the exact ground-state local one- 
body density po(r) = ('F 0 |a+a r |'I , o)- It is also obvious 
that the above argumentation can be repeated mutatis 
mutandis for a functional of a one-body non-local density 
p(r,r') = (Wla^/Or|’P), which is the formulation we are 
concerned with below. 

Certainly, such an idealistic derivation of the exact 
EDF would defy its purpose: had we been able to perform 
the exact variational calculations for all one-body poten¬ 
tials —U(r), we would have probably not need DFT at 
all. The strength and beauty of DFT is somewhere else: 
general considerations about the exact DFT provide us 
with a motivation to search for a suitable and physically 
justified modelisation E\p\ of the exact EDF E\p\. In 
this way, one is only left with an easy task of performing 
the second-stage variation over the density p{r). Unfor¬ 
tunately, the rigorous link between the exact many-body 
Hamiltonian H and the model EDF E [p] is then lost. In 
this paper, I propose a method of recovering it. 

The goal is thus not to derive the exact EDF E [p] , but 
to provide an ab initio derivation valid within a certain 
class of model EDFs E \p\. Such models should be spe¬ 
cific to a given range of energies or distances, at which 
low-energy description of ground states of given physical 
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systems is relevant. 

The class I am going to employ is motivated by 60- 
odd years of modelling EDFs in nuclei j3l[, and can be 
formulated as 


= ( 2 ) 

»=i 


where C l are coupling constants and Vi [p] are the 
Hartree-Fock (or first-order many-body-perturbation- 
theory) averages of certain two-body, three-body, etc., 
operators Vj. At early stages of developing the nuclear 
EDFs, these operators were called interactions, but in 
fact, their sole role was to generate specific terms in the 
EDF, so here I call them EDF generators. 


For the construction presented below, it is essential 
that the model EDFs (12) are built in terms of true oper¬ 
ators acting in the many-body space, because one must 
be able to use them not only for defining the EDFs, but 
also within the true ab initio many-body context. On 
the one hand, some constructs typical in nuclear EDFs, 
like the explicit density-dependent terms .3lj, are thus 
excluded. On the other hand, functionals based on EDF 
generators seem to be the only ones that allow for us¬ 
ing EDFs in the multi-reference context, see, e.g., recent 
Ref. [32|, and, therefore, constructions based on EDF 
generators are very much called for. We note here that 
the proposed scheme would also work for EDFs generated 
by operators depending on additional parameters, so the 
specific linear dependence on the coupling constants al¬ 
though convenient, is not really essential. 


Before considering specific EDF generators Vi that 
were used and/or proposed in nuclear physics, let us dis¬ 
cuss the main consequences of using the model EDF in 
the form of Eq. ©• First of all, one should keep in mind 
that the EDF is always meant to be minimized with re¬ 
spect to the density, and thus its detailed form beyond 
the minimum is not essential. By the same token, there is 
always a one-to-one correspondence between the coupling 
constants of the functional C l and densities that mini¬ 
mize it. Therefore, the manifold of meaningful ground- 
state densities M[p\ is not really infinite dimensional, but 
it can be parametrized by the coupling constants C\ and 
eventually by conserved quantum numbers, so it has a fi¬ 
nite number of dimensions. Conversely, the model EDF 
© does not have to properly describe the exact ener¬ 
gies of states having all possible densities, but only those 
that have densities on this restricted finite-dimensional 
manifold M[p\. 


This important observation has far reaching conse¬ 
quences. Indeed, instead of probing the system with 
all possible one-body potentials — U(r ) of an arbitrary 
shape, as in Eq. ©, it is enough to probe it within the 
finite set of the EDF generators — Vj, that is, to solve the 


constrained variational equation, 

m 

6E' = 6($|F-^A- 7 'y 7 -|'F) =0, (3) 

t=i 

for a suitable set of values of a finite number of Lagrange 
multipliers A*, which is perfectly manageable a task. In 
Eq. ©, there appear the same EDF generators, which in 
Eq. © were used to define the model EDF in the first 
place. This is perfectly logical: to meaningfully include 
a term in the model EDF we must first test its proper¬ 
ties in the real world of the ab initio phase space and 
Hamiltonian. 

Solution of Eq. © gives us the exact ground-state en¬ 
ergies E( A- 7 ) and one-body non-local densities p X j (r i, r 2 ), 
both as functions (not functionals!) of the Lagrange mul¬ 
tipliers A 7 . Of course, now the dependence of densities on 
Lagrange multipliers cannot be inverted, however, this is 
not at all necessary. It is enough to ensure that, on the 
manifold generated by the Lagrange multipliers A- 7 , the 
model EDF © best reproduces the exact energies, that 
is, it is enough to adjust the EDF coupling constants C l 
so as to have, 

m 

E{X) = Y J C i V i [p xi }- ( 4 ) 

i=1 


The adjustment is performed for a finite set of values 
of the finite set of Lagrange multipliers, so Eq. © con¬ 
stitutes, in fact, a basic standard linear-regression prob¬ 
lem. After the adjustment, one obtains a true ab initio- 
equivalent EDF. 

The ab initio derivation of the model EDFs, proposed 
in this work, may become a basis for future studies that 
can bridge the ab initio methods with those related to de¬ 
riving and improving the phenomenological EDFs. The 
proposed research program will probably take some time, 
especially in view of the fact that present-day successful 
ab initio implementations are at the forefront of what 
is currently possible within the high performance com¬ 
puting. Therefore, in this work I only present a simple 
proof-of-principle application of the proposed scheme to 
a task of relating one class of the EDF to another. 

To this end, I used the EDF generators corresponding 
to the central and tensor parts of the nuclear Skyrme 
interaction 3l|, [33-36J, which is composed of eight terms 


(m=8), that is, 


M \ 


( 4 

-2 

0 

0 

0 

0 

0 0 

\ 

(T 0 \ 

K 


0 

-6 

0 

0 

0 

0 

0 0 


rp<7 

V 0 Ap 


0 

0 

-8 

4 

8 

-4 

00 

1 

00 


Ti 

vf p 

2 

0 

0 

0 

12 

-16 

20 

-24 -8 


fa 

V T 

v o 

~ 3 

0 

0 

2 

-1 

6 

-3 

4 0 


T 2 

V{ 


0 

0 

0 

-3 

-12 

15 

-12 0 


T% 

V 0 J 


0 

0 

0 

0 

0 

0 

-6 -6 


T e 

W / 



0 

0 

0 

0 

0 

18 -6 

/ 

UJ 


( 5 ) 














3 


where Tq = 5, T\ = ^(fc'* 2 + k 2 )S , T 2 = (k 1 * • fc)<5, and 
(5 is a two-body local zero-range potential, 5 = 5(ri — 
i , '2)S(i'i—r' 1 )5(r2—r! 2 ). The standard relative-momentum 
operators are defined as k = (Vi — V 2 V 2 i and k' = 
(Vi - V' 2 )/2 i, T? = i\P G for Pfj = ±(1 + <xi • a 2 ), and 
T e = \{k'* • S • k 1 * + k ■ 5 ■ k) and T 0 = k'* • S ■ k for 

S ab _ 3 (—a —b i _6 ~a\ v _ _ 

— 2H^2 ' ) Oab& 1 ' • 

Numerical coefficients appearing in Eq. © were chosen 

in such a way that each of the eight EDF generat ors give s 

(in spherical nuclei) one specific term of the EDF |36H38j . 

namely, 

V t p \fi] = pt{r)p t (r) , V t Ap [p\ = pt(r)Ap t (r), 

V t T [p] = Pt{r)Tt(r) , V t J [p] = J t (r) • J t (r), 


where index t refers to the isoscalar (f=0) or isovector 
(£=1) densities, Ap t (r) stands for the Laplacian of the 
density, and r t (r) = [(V • V')/9 t (r, r')] r _ r , and J“ b (r) = 
A [(V a — V(Js b (r, r')] are the standard quasilocal 
kinetic and spin-current densities, respectively. 

For the proof-of-principle application presented in this 
work, instead of the average value of the true many- 
body Hamiltonian, I used the Gogny EDF j3f| in the 
D1S parametrization , that is, ->• Ecogny [p\- 

In this way, I aimed at obtaining a Gogny-equivalent 
quasilocal Skyrme EDF. That both types of EDFs can 
be linked to one another is already known from findings 
of Refs. iHIH , where this fact was demonstrated within 
the DME and effective theory; here I aim at testing this 
equivalence in terms of the ab initio -like methodology 
proposed in this work. Since both functionals contain 
terms generated by the zero-range spin-orbit and density- 
dependent operators, these were left untouched, and in 
the left-hand side of Eq. only the part of the Gogny 
EDF generated by the finite-range potentials was used. 

Numerical results presented below were obtained using 
the code hfodd (v2.75c) [43|> which is the only existing 
code capable of treating the Gogny and Skyrme function¬ 
als simultaneously and within the same numerical infras¬ 
tructure, see the Supplemental Material [H| for details. 
Calculations were performed for eight doubly magic nu¬ 
clei, 16 0, 40 ’ 48 Ca, 56 ’ 78 Ni, 10 °. 1 32g n , and 208 p k For each 


nucleus, I used either of the eight Lagrange multipliers, 
Aj , X Ap , X[, or A/ (for t=0,l) equal to one of the 21 in¬ 
teger values between —10 and +10MeV fm", where n =3 
for X p and n =5 for the other ones. Altogether, this gave 
me 1344 values of the Gogny energies E( A J ), to which the 
eight coupling constants of the Skyrme EDF, C p , C t Ap , 
C t T , and C f J (for t=0,l) were adjusted in Eq. ©■ 

The coupling constants SISd obtained by such an ad¬ 
justment are shown in TablelU Their standard uncertain¬ 
ties were obtained within the standard regression analysis 
presented, e.g., in Refs. jiH, J46[. The adjusted coupling 
constants reproduced the Gogny energies in Eq. © with 
a very high accuracy: the relative rms deviation between 
left- and right-hand sides of Eq. © is only 0.014%. This 


TABLE I: The Skyrme EDF SISd and SISe coupling con¬ 
stants obtained in this work as the ab initio-equivalent Gogny 

EDF D1S [IJ. 

SlSd STSi 

i=0 t=l t= 0 t=l 

C t p (MeVfm 3 ) -603.82(22) 484(4) -605.41(16) 509(3) 

C t Ap (MeVfm 5 ) -73.25(14) 48(3) -74.82(12) 41(2) 

Cl (MeVfm 5 ) 77.95(23) -78(3) 79.73(16) -98(2) 

Cj (MeVfm 5 ) 24.9(1.2) 71(3)_0 0 



FIG. 1: (Color online) Gogny energies [lines, left-hand side of 
Eq. Q] compared with the EDF estimates [symbols, right- 
hand side of Eq. (©] obtained for the Skyrme EDF SISd 
coupling constants given in Table [I] Calculations were per¬ 
formed in 208 Pb in function of the eight Lagrange multipliers 
A p , A^ p , A l, and A/ for £=0,1. The inset shows residuals of 
the adjustment in per cent. 


is also illustrated in Fig. [I] where differences between 
symbols and lines cannot be seen at all, while the inset 
shows that the relative residuals of the adjustment do 
not exceed 0.050%. Analogous plots for other nuclei are 
collected in the Supplemental Material [ddj. 

Table HD compares the ground-state energies Eg calcu¬ 
lated using the original Gogny EDF D1S with energies 
E obtained by the minimization of the Gogny-equivalent 
Skyrme EDF SISd. Propagated uncertainties A E of E 
were calculated using the covariance matrix related to the 
adjustment of coupling constants HHEsI, see Supple¬ 
mental Material [3]. We see that the Skyrme EDF SISd 
again reproduces the Gogny-EDF results with a very high 
accuracy: the relative rms deviations between these two 
functionals is only 0.28%. This is much better than the 
accuracy obtained within the DME [47j, see the cornpar- 
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TABLE II: Gogny EDF D1S ground-state energies Eg (b) 
of eight doubly magic nuclei (a) compared to energies E (c) 
calculated using the Skyrme EDF SISd, Table Q] and shown 
together with their propagated uncertainties A E. Column (d) 
shows the residuals 5E = E—Eg and columns (e) and (f) give 
ratios of residuals with respect to energies E and propagated 
uncertainties A E, respectively. All energies are in MeV. 


(a) 

Eg 

(b) 

E 

(c) 

5E 5E/\E\ 
(d) (e) 

SE/AE 

(f) 

-TOq 

-129.626 

-129.56(4) 

0.07 0.05% 

2 

40 Ca 

-344.663 

-346.01(6) 

-1.35 -0.39% 

-23 

48 Ca 

-416.829 

-418.10(6) 

-1.27 -0.30% 

-20 

56 Ni 

-483.820 

-485.31(6) 

-1.49 -0.31% 

-27 

78 Ni 

-640.598 

-642.37(6) 

-1.78 -0.28% 

-28 

100 Sn 

-830.896 

-833.19(6) 

-2.29 -0.28% 

-39 

132 Sn 

-1103.246 

-1106.51(8) 

-3.26 -0.30% 

-42 

2° 8 p b 

-1638.330 

-1640.96(8) 

-2.63 -0.16% 

-32 

rms 

n.a. 

n.a. 

1.99 0.28% 

29 


ison presented in the Supplemental Material [44], One 
can conclude that the simple eight-dimensional Gogny- 
equivalent Skyrme EDF with ab initio- derived coupling 
constants of Eq. (0 very well describes the full Gogny 
energies. We note here that the coupling constants were 
adjusted only to energies. In the Supplemental Mate¬ 
rial 44j, I also show the analogous very good agreement 
obtained for the proton rms radii. This points to a well- 
built EDF, which along with the total energies properly 
describes one-body observables. 

On the absolute scale, the corresponding rms deviation 
of energies is 1.99 MeV, which is 29 times higher than the 
rms average of the propagated uncertainties, see Table HT1 
It means that the differences between the Gogny results 
and Gogny-equivalent Skyrme-EDF results are still sig¬ 
nificantly larger than the uncertainties of the adjustment, 
which gives a clear signal for missing terms in the eight¬ 
dimensional model EDF of Eq. ©• 

The method also allows for testing the impact of re¬ 
moving terms from the model EDF. For example, set¬ 
ting the two spin-current coupling constants to zero, 
c 0 j = a = 0, one obtains a six-dimensional model 
that gives another Gogny-equivalent Skyrme EDF SISe, 
with coupling constants shown in Tables [U Such model 
is only marginally worse, with the absolute and relative 
rms deviations of energies now increased to 2.34 MeV and 
0.40%, respectively, see the Supplemental Material [441 
for detailed results. 

In conclusion, I proposed a novel method of obtaining 
ab initio-e quivalent model EDFs. The main idea is in 
replacing the standard-DFT use of an external one-body 
potential by the use of two-body, three-body, etc., EDF 
generators. This probes the reaction of the system with 
respect to the same operators that are used to construct 
the model EDFs. The new method amounts to perform¬ 
ing ab initio calculations with simple constraints on a 


finite set of well-defined operators added to the many- 
body Hamiltonian, and thus is perfectly manageable. 

The method is also able to give us a quantitative in¬ 
formation on whether a given model EDF is adequate 
for the proper description of the physical system being 
studied. Indeed, if the adjustment of the coupling con¬ 
stants fails to be accurate enough, we obtain a clear sig¬ 
nal that the set of proposed EDF generators is inade¬ 
quate. Then, another or extended model EDF should 
be tried. In this way, through ab initio derivations in 
nuclei, one may be able to evaluate relative merits of 
using the zero-range higher-order [48|, finite-range reg¬ 
ularized three-bodv or four-body [50], or gradient- 
dependent three-body [5JJ] pseudopotentials, which are 
currently being developed and implemented in nuclear 
EDF approaches. 


It is very important that the proposed method is based 
on studying specific finite systems and does not rely on 
assumptions valid only in infinite or semi-infinite sys¬ 
tems. The proposed ab initio derivations can be per¬ 
formed in few systems, for which the ab initio calcula¬ 
tions are possible, e.g., in light closed-shell nuclei, and 
then the derived EDFs can be applied to more compli¬ 
cated, open-shell or heavy nuclei, so as to test the over¬ 
all predictive power of the method. It is also essential 
that the ab initio-equivalent EDFs are specific to par¬ 
ticular physical systems, and when applied to systems 
at different energies or densities will yield parametrically 
energy- or density-dependent running coupling constants. 
Needless to say that the proposed method can easily be 
extended to deriving time-odd or pairing terms of the 
EDFs. Another fascinating extension would be to use the 
same method not only to match ener gies , like in Eq. 
but also ab initio- derived kernels 52), [53J. 
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